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At coupling strengths A = 1/2, 1, or 2, the Calogero-Sutherland model (CSM) is related to Brow- 
nian motion in a Wigner-Dyson random matrix ensemble with orthogonal, unitary, or symplectic 
symmetry. Using this relation in conjunction with superanalytic techniques developed in mesoscopic 
conductor physics, we derive an exact integral representation for the CSM two-particle Green's func- 
tion in the thermodynamic limit. Simple closed expressions for the single-particle Green's functions 
are extracted by separation of points. For the advanced part, where a particle is added to the 
ground state and later removed, a sum of two contributions is found: the expected one with just one 
particle excitation present, plus an extra term arising from fractionalization of the single particle 
into a number of elementary particle and hole excitations. 
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I. INTRODUCTION 



While the study of integrable models of interacting 
particle systems in one dimension has a long history, go- 
ing back to Bethe j^] , it has proved difficult to extract ex- 
plicit expressions for their correlation functions from the 
solutions. Recently, however, the correlation functions of 
one class of such models, those with inverse-square inter- 
actions, have started to become available. These models 
have a special simplicity as compared to the integrable 
models of the Bcthc type, and combine some of the spe- 
cial properties of an ideal gas with non-trivial statistics. 

The original model of this type is the Calogero- 
Sutherland model §,| (CSM) of a gas of spinless non- 
relativistic particles on a ring, with a scale-invariant in- 
teraction proportional to the inverse-square chord dis- 
tance between particles. Different members of the CSM 
family are characterized by a dimensionless parameter 
A, where the wavefunction vanishes as {xi — Xj)^ as two 
particles approach each other. At the coupling A = 1 the 
interaction vanishes, and the model reduces to a gas of 
free spinless fermions. Early progress in obtaining cor- 
relation functions of the model came from Dyson's work 
on correlations of the eigenvalues of random matri- 
ces: at the special couplings A = 1/2, 1, and 2, where the 
CSM is related to the orthogonal, unitary, and symplectic 
random matrix ensembles. Dyson's work led to explicit 
expressions Q for the static density-density correlations 
at A = 1/2 and 2, and for the density matrix at A = 2. 
For many years, these were the only explicitly-obtained 
correlations of integrable models. 

The recent breakthrough derives from the work of Si- 
mons, Lee and Al'tshuler [| (SLA), and again has its 
origin in the study of random matrix problems. Using 
the functional integral formalism introduced by Efetov 



Q, SLA obtained expressions that they initially conjec- 
tured were the dynamical density-density correlations at 
A = 1/2 and 2, since the analogous quantity at A = 1 
was indeed that quantity for the free Fermi gas, and the 
equal-time limit of their result coincided with the results 
known from Dyson's work. Their remarkable conjecture 
was later verified. A variant of the SLA methodology was 
used by us (HZ) to obtain the retarded part of the 
bosonic single-particle Green's function at A = 2, which 
at equal times reduces to the previously-known expres- 
sion for the density matrix. When rewritten in terms of 
form factors^ the SLA results become especially simple, 
and led to a conjecture by one of us Q of their general- 
ization to all rational coupling parameters A. This con- 
jecture was recently proved using Jack polynomial tech- 
niques Q (and was independently made at integral values 
of A, using group-theoretical arguments |]lO|). Jack poly- 
nomial methods were introduced in this context by For- 
rester PI , who recently jl^ used them to prove a similar 
conjecture [|l^ extending the HZ result to all integer A. 
Ha [|j has also extended this to arbitrary rational A. 

While Jack polynomial methods are not restricted to 
particular values of A, they require further development. 
In the meantime, the random matrix techniques provide 
powerful computational methods at couplings A = 1/2 
and 2. In this paper, we extend our earlier result |j] 
by also giving the advanced part of the bosonic single- 
particle Green's function at A = 2, and extending the 
calculation to the other coupling A = 1/2. 

The paper is organized as follows. Sect. || reviews 
the physics of the CSM, with particular emphasis given 
to the identification of its elementary particle and hole 
excitations. Anyonic particle creation and annihilation 
operators are defined in Sect. III. Our results are stated 
in Sect. |l^. The two-particle Green's function defined in 
Sect. is rewritten as a harmonic oscillator correlation 
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function in S ect. VI. The bulk of the paper is contained 
in Sect. VII , where the supersymmetric calculation of 
Green's functions is described. Many technical details 
have been relegated to three appendices. 



II. THE CALOGERO-SUTHERLAND MODEL 

The Calogero-Sutherland model is a model of 

a gas of non-relativistic impenetrable particles with 
inverse-square interactions. It is convenient to impose 
periodic boundary conditions on a ring of length L, al- 
though the simple results for the correlation functions 
will only be given in the thermodynamic limit. The 
Hamiltonian is 



H, 



cs 



^ 2m 



E 



d{xi - XjY 



(1) 



where d(x) = (L/tt) sm{Trx/ L) is the chord distance be- 
tween particles on the ring. The model is fully defined 
by the boundary condition on the wavefunction when two 
particles approach each other: 



1*1 oc \x, 
where A > 0, and 



0, 



9 = • 

TO 

Thus A ^ corresponds to the free Bose gas, and A ^ 1 
to the free (spinless) Fermi gas. (For —1/2 < g < 0, 
there are two different choices of A for each g.) 

The energy spectra |^ are given by the Bethe-Ansatz 
like equations, and can be viewed as adiabatic deforma- 
tions of the free-electron spectrum for A = 1: for i — 
1,...,7V 



kiL = 27r/i + 7r(A — 1) sign(fci — kj), 



(2) 



where the {li} are distinct integers (or for more general 
boundary conditions, all h — Ij are integral). Energy and 
momentum of the eigenstate are then given by 



(3) 



In the ground state, the /i's are consecutive integers, 
and the fc^'s are equally spaced, with spacing 2ttX/L. 
If an excited state is made by creating a "hole" in this 
"pseudo- Fermi sea", by removing one of the /^'s, the gap 
between the two consecutive ki on either side of the hole 
is 27r(A-|-l)/L = (27rA/L)(l -h A'^). There is a very sim- 
ple interpretation of this: the dressed charge of the hole is 
— A^^, as the extra width of the pseudo-Fermi sea associ- 
ated with the presence of the hole is A""'^ of the width per 
particle. (Here the bare particles are taken to have unit 



charge). The other kind of excitation is to add a particle 
in a momentum state outside the pseudo-Fermi sea. In 
this case, the charge is -1-1, the bare charge. In the free 
fermion limit, A = 1, these identifications coincide with 
the usual ones for the Fermi gas. 

Besides carrying fractional charge, the excitations have 
a natural interpretation as particles carrying fractional 
statistics (see e.g., |p^ ). This is perhaps easiest to see 
from the wavefunctions for the periodic model, which 
have the form Hi 



(4) 



where Zi = e^'^"*/^, 0({zj}) is a symmetric polynomial 
that is not divisible by Zi (called a Jack polynomial |^ ) , 
J is any integer, and a is fixed (modulo an integer) by the 
choice of generalized periodic boundary condition (the 
phase change of 5* when Xi Xi + L) . In the older liter- 
ature 1^, it was usual to replace (zi — zj)^ by \zi — Zj\^ 



or Zi 



and call the bare particles bosons 
or fermions, but in the light of the striking resemblance 
of the CSM wavefunction to the Laughlin wavefunction 
iQ, and the wavefunction for "anyons", the form (^) 
seems more appropriate. This strongly suggests that the 
bare particles of the CSM be identified as anyons, with 
"statistical parameter" Q = nX. Polychronakos |l^] has 
also given arguments for this based on scattering phase 
shifts. 

The model for anyons is that of particles carrying both 
charge and magnetic flux jl^]. We thus can identify 
the bare CSM particle as effectively carrying "electric" 
charge Q and "magnetic" flux $, so the statistical angle 
6 is Q^/2fi, where 

particle: Q = e, $ = A/i/e, 9 = ttA. 

What about the hole excitation? A "coherent" state 
with a single hole has the wavefunction 



where ^'o is the ground state wavefunction. Here Z (if 
unimodular) parameterizes the position of the hole; ex- 
pansion of this wavefunction in powers of Z gives its ex- 
pansion in states of definite momentum which are the 
one-hole eigenstates of the Hamiltonian. All this is es- 
sentially the same as in the Laughlin states and the frac- 
tional quantum Hall effect (FQHE) Q. Taking over the 
charge-statistics results from the FQHE ||l^ , leads to the 
identification of the charge Qh, flux <i>h, and statistical 
angle 8h of the hole: 



hole: Qh 



-X- 



-h/e, Qh ~ 



Again, these identifications reduce to the trivially- 
expected ones when A = 1, but are otherwise highly non- 
trivial. 
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The CSM has Gahlean invariance, and it is straight- 
forward to compute the effective mass of the holes from 
the equations (The inertial mass of the particle 

excitation is unchanged from the bare value). The hole 
has inertial mass mh — —mX~^, so that, as expected in 
a Galilean- invariant system, Qh/^^-h = Q/m. This is fur- 
ther corroboration of the identification of the fractional 
charge of the hole. Since the particle and hole excita- 
tions have different masses, the appropriate dynamical 
variable with which to describe them is their velocity. 
Holes have velocities in the range —Vg < v < Vs, where 
Vs = TrXhpo/m is the speed of sound (long wavelength 
density fluctuations) in the model, and po — N/ L is the 
mean density. Left-moving low-energy holes w « and 
right-moving low-energy holes v ~ Vs are continuously 
connected to each other. On the other hand, the particle 
excitations have \v\ > Vs, and the left- and right-moving 
states are disjoint. 

If an eigenstate is characterized by a set of particle ex- 
citations with velocities {vi} and a set of hole excitations 
with velocities {vi}, momentum and energy (relative to 
the ground state) are given by 

P = '^mvi + ^TOhWi, 

i i 

E^Y. - + E l"^h(^;? - vl). (5) 

i i 

Let us now consider the implication of these identifica- 
tions for the form factors: what types of excitation are 
produced when the density operator 

p{x) = Y^S{x - Xi) 

i 

acts on the ground state? The density operator carries 
no charge, and is a bosonic operator. For A = 1, we 
know the answer: the action of p(x) on the ground state 
produces exactly one hole and one particle, with a total 
charge Q + Qh — Q, and a total flux $ -I- $h = 0. 

The recent calculations of Simons, Lee, and Altshuler 
when suitably interpreted, give the non-trivial an- 
swer for A = 2 and A — 1/2. For A = 2, one particle and 
two holes are produced, while for A = 1/2, two particles 
and one hole are produced. This count corresponds to 
the number of free velocity parameters in the form fac- 
tors implied by the SLA results. As expected, these are 
consistent with a total charge of zero and a total flux 
of zero. The generalization of this is that at rational 
coupling A = p/q, the action of the local density opera- 
tor on the CSM ground state will produce excited states 
with exactly q particles and p holes. A further restric- 
tion is found: in the A = 1/2 CSM, the results ||] show 
that both the particles have velocities in the same di- 
rection (momenta on the same side of momentum space 
with respect to the pseudo- Fermi sea). The generaliza- 
tion of this is the result that all the q particle excitations 
have velocities in the same direction (all left-movers or 



all right-movers). Thus the simplest possible excitations 
consistent with charge neutrality account for the com- 
plete spectral weight of the density-density correlation 
function. This is a generalization of the free Fermi gas 
property of the A = 1 model, where just one particle and 
one hole are produced by the action of the local den- 
sity operator. In a more general interacting model, there 
would be contributions from two particle-hole pairs, three 
particle-hole pairs, and so on. The remarkable structure 
of the non-trivial solutions found at A = 1/2 and 2 is 
that the only intermediate states that contribute are the 
minimally-excited ones, where the smallest number of ex- 
citations consistent with charge neutrality are made. In 
this sense, the CSM behaves like a generalized ideal gas 
model. 

In the case of the single-particle Green's function, there 
arc two terms to consider: the advanced term, where a 
particle is added to the system, and removed at a later 
time, and the retarded term where it is removed and 
later restored. In these cases, the intermediate states 
must again have an excitation content with the same to- 
tal charge and flux as the operator that initially acts. 
Thus in the advanced case, the total charge is e, and the 
flux is Xh/e; the negatives of these apply in the retarded 
case. At coupling A — p/q, we find the expected contribu- 
tion with just one particle excitation, plus a second term 
with one particle in the forward direction, q particles in 
the backwards direction, and p holes. In the retarded 
case, the general result found by Ha is q — 1 particles (in 
the same direction) and p holes. 



III. ANYON CREATION AND ANNIHILATION 
OPERATORS 

We have reviewed the arguments in favor of the inter- 
pretation of the CSM particles as anyons when A is not 
an integer. This is the interpretation we adopt. 

To calculate the anyonic CSM Green's functions, we 
shall need to define the action of the particle creation 
and annihilation operators on A'^-anyon states. Recall 
first how the particle annihilation operator ip^ acts on a 
state of N bosons (s = +1) or fermions (s — —1): 

N 

{{Xi}\tpx\'^ n) = y^^s'^Nixi,...,Xi-i,X,Xi+i,...,XN). 

1=1 

Now, the Hilbert space of A^-anyon CSM wavefunctions 
is spanned by states of the form 

{{x^}\^n) ^ SNi{x^})Y[dix, ~ x,)\ (6) 

i<j 

where SN{{xi}) is a symmetric function. When A is 
an even (odd) integer, these wavefunctions are bosonic 
(fermionic), and in these cases insertion into the previ- 
ous equation gives (up to normalization): 
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i<j 

^'N,li{^i}) = 5Ar+i(x, {a;^}) ]^d(a; - Xi)^. 



(7) 



Note that the function S*^ i is a symmetric function of 
its N arguments, not only for integral A but for any A. 
Therefore, it makes sense to define the annihilation op- 
erator V'a; by (|7|) at any A > 0, and this is what we do. 
The operator so defined properly maps the (A'^-l-l)-anyon 
Hilbert space onto the 7V-anyon Hilbert space. Similarly, 
by starting from the formula 

N+l 

{{Xi}\-ipl\^ n) = ^ S''5{x - Xi)^Ni-;Xi-i,Xi+i, ...) 
i=l 

for bosons or fermions, we obtain 
{{x^Ml^M-i) = S^+l{{x^.})Y[d{x, - x,)\ 



5jv-i({a;-,, j 7^ i})5(x - Xj) 



Again, the last equations make sense, and give a valid 
definition of the anyon creation operator i/jI, for any A. 



IV. STATEMENT OF RESULT 

Our results are as follows. Let |0) denote the A^- 
particle ground state. In the thermodynamic limit N ~* 
oo, po = N/L fixed, the advanced part of the single- 
particle Green's function {t > t'), 

Gp(x-x';t-t')-(0|V^Jt)^l,(t')|0), 

separates into two terms: Gp ~ Gp^' -I- Gp^^ . To present 
the result for the cases A — 1/2,1,2 simultaneously, we 
set A ~ p/q. The first term has the elementary form 



Up (x,t) 



V + Vs 



cos(fca;)e-'"*du 



where k — mv/h, uj — m{v'^—Vg)/2h, and Vg — irXfipf^/m. 
The second contribution arises from fractionalization of 
the single particle into g-|-l particle and p hole excitations 
with velocities V := {vi, Wq+i} and V := {vi, Vp\: 



(8) G^^\x-t)= \F{q+l,p,\\V,V)? cos{kx)e-"^^dVdV. 



Here dV = 11?=! ^^i^ = n?=i dvj, and the integra- 



tion is restricted to the domain 

OO < Vl < ... < Vq < —Vg < Vl < ... < Vp < +Vs < Wq+l < -|-00. 



Wave number k = P/Ti and frequency uj = E/Ti have been given in 
consists of two parts: F = F^^'> x F^^) . The first one. 



The "form factor" \F\'^ = \{V,V\'ipl,\Q)\ 



|F(°)(np,nh,A|y,y)|' = 



(9) 



arises in our calculation as a jacobian or measure and in 
this sense is purely statistical, doing no more than mea- 
suring the phase space available to the elementary excita- 
tions, is completely determined by the root system 
of an underlying Lie superalgebra and is therefore uni- 
versal. (More precisely speaking, if rip and rih are the 
numbers of particles and holes an excitation is composed 
of, the form factor of that excitation will contain a factor 
F("'(np,nh, ...).) Note also the duality symmetry 

F(°)(np,nh,A|-l/,F) =F(")(nh,np,A-i|-^,F). (10) 

The other factor is specific to the advanced part of the 
single-particle Green's function and is given by 




The normalization constants are G2 = y^m/8Trh and 
Gi/2 = y/m/4:TT^h. 

Our calculation of F^^ for A = 1/2,1,2 is suggestive 
of generalization and leads us to propose the following 
conjecture for integral A = p > 1: 



F 



(1) 



Cx I -1- 

dv 



{v — Vi)^ 



(11) 



Verification of this conjecture is posed as a challenge to 
the theory of Jack symmetric functions. 



V. TWO-PARTICLE GREEN'S FUNCTION 



The objective of this paper is the calculation of the ad- 
vanced part of the single-particle Green's function. For 
reasons of generality and convenience of presentation, we 
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will set up our computational machinery for the two- 
particle (or particle-hole) Green's function (t > t') 



,(t-t') = (0|^J(t)V'.W^l(i')'A.(i')|0). 



For X — y and w = z, this Green's function reduces to 
the dynamical density correlation function calculated in 
When the pair x, w is taken to be remote from the 
pair y, z, we expect G to separate: 

Gxywz{i) ^ Gp{x - w] t)Gh{y - z; t) 

where Gp and GY,{x-y]t-t') = {'^\'4>lit)iljy{t')\{)) are the 
retarded and advanced parts of single-particle Green's 
function. 

Our first task is to describe the action of the particle- 
hole creation operator V'lV'b on the A'^-particle ground 
state |0). By combining the definitions (|^, ||) we see that 
this action is simply multiplication with a function of the 
coordinates X :— {xi, ...,X]y}: 



(X|VlA|0)=O,,(X)(X|0), 



N 



Oab{X)=J2S{a-x 



d\b~xk) 

^^n.d^-.,)- 



)n 



(12) 



As it stands, the "particle-hole function" Oab{X) is not 
in a form suitable for calculation by the technique used 
in the present paper, which is why we make the following 
modification. Consider for e > the function : R — > R 
defined by (3? is the real part) 

Seix) = 3? (2ie)^(27ri)-i(.T + is)-^{x - ie)-\ 

Since /(O) = lim^^o Jn fix)6e{x)dx for any smooth 
function /, is a valid regularization of Dirac's delta 
distribution centered at zero. We make 6^ periodic by 
replacing distance by chord distance: 



d^{x + i£)d{x ~ ie) 
Inserting this into the expression for Oab{X), we obtain 

o.bix) = ^jm^. 



d-^(6 — a + ie) 



N 

E 



N 

n 



d^(6 - Xk + is) 
^ d(a — Xj — is) d^(a — Xk + is) 

1 k—l 



(13) 



which differs from (^^ by terms that become negligible 
in the limit e ^ 0. Introducing 

K{X,X';t) = {0\X){X\e-''"'^^/''\X'){X'\0) 

we can now write the two-particle Green's function in the 
form 

Gxy^,{t)= JdX JdX' Oxy{X)K{X,X':t)0^,{X'). 

Here, dX — Y[ dxi and the integration range is that sec- 
tor of [0, L)^ where Xi < Xi+i for i = 1, TV — 1. 



VI. MAPPING ON A SYSTEM OF HARMONIC 
OSCILLATORS 



It has long been known ||T^ that the CSM at coupling 
strengths A = 1/2, 1, and 2 is equivalent to free-particle 
motion on a symmetric space U(A^)/0(iV), U(iV) and 
U(2A^)/Sp(2iV), respectively. More precisely speaking, 
for A = 1/2, 1, or 2 there exists a similarity transforma- 
tion taking i?cs into the radial part of the kinetic energy 
operator on the corresponding symmetric space. We will 
now review this transformation for the case A = 2. 

Let C = l_N be the symplectic unit acting on the 

2A^-dimensional linear space (g) C^. With g running 
through the group U(2iV) we consider the set, Mn, oi 
matrices S of the form 

Sig) = gCg^C-\ 

Such matrices are unitary and obey the self-duality con- 
straint S = CS"^C~^. They do not form a group. Rather, 
since S{gk) = S{g) for any element k of the symplec- 
tic subgroup Sp(2A^) C U(2A^) defined by the condi- 
tion kCk'^ = C, Mn isomorphic to the coset space 
U(27V)/Sp(27V). U(2A^) acts as a transformation group 
on Mn by S{g) i-> S{hg) {h e \J{2N)). When endowed 
with its natural U(2A^)-invariant metric — ^Trd5d5^ 
and probability measure dfi{S), A4n is known in random 
matrix theory as Dyson's Circular Symplectic Ensemble 
(CSE) |,|l|. 

Dyson, in his classic paper proved the following 
two results, (i) Any element S Q Mn has a "polar" 
decomposition S = ke^^k~^ where k g Sp(27V), and 
9 = dia.g{6i, ...,9n) <Xi I2 is diagonal. (We refer to e'^ 
and k as the "radial" and "angular" parts of S. Note 
that every eigenvalue of S occurs with multiplicity 2.) 
(ii) Polar decomposition takes the invariant integral on 
Mn into 



fiS)dfiiS) 



[0,27r)« \ JSp(2Af) 



where d6 — Ht^i'^^i' ^nd dk is a suitably normalized 
Haar measure for Sp(2Ai"). The jacobian J for the trans- 
formation to polar coordinates is 

j{e) = l[sm'[ie,-e,)/2]. 

Now, let C be the Laplacian (or minus the operator for 
kinetic energy) of the Riemannian manifold Mn- If ^Cr 



f{ke''^k^^)dk J{e)d6 
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is the radial part of C, we have £/ = C^f for any radial 
function /. Using 



N 



TrdS-d^t = + Ti[k-^dk,e''^f 



in conjunction with the standard coordinate expression 
for the Laplacian on a Riemannian space, we easily get 



i—l 

The differential operator is self-adjoint w.r.t. the ra- 
dial measure J{6)d0. It is made self-adjoint w.r.t. the flat 
measure dO by carrying out the similarity transformation 
A- Ji/2£^j-i/2^ ^j^-jj tl^g jjglp ^-jjg identity 



WijWjk + WjkWki + WkiWij = 1 

cot[(6'i — 6'j)/2], and by making the identifica- 



for Wi 

tion 6i = iiTXi/L, we obtain 

This is the desired relation between the CSM at A = 2 
and free-particle radial motion on the symmetric space 
Mn ^ U(27V)/Sp(27V). 

For A = 2, the function Oah{X) of Eq. (|l|) can be 
viewed as a radial function 



27F^ l-K 



d2(6 — a + ie) 



-Tr T- 

2 d(a-^ 



Det 



d(6- 2|jln5 + z£) 



InS'-ie) d(a- 2^1nS' + ie) 



on A^AT. Hence, the particle-hole Green's function of the 
CSM can be calculated as a dynamical correlation func- 
tion for free-particle motion on M-fq. Note that compu- 
tation of the latter is not entirely trivial, since A^at is a 
curved space. 

A further simplification arises in the thermodynamic 
limit — > oo. In this limit, Dyson's circular ensem- 
bles are known to become locally equivalent to the corre- 
sponding Gaussian ensembles [ppS|]. For the dynamical 
correlation functions, the equivalence can be understood 
by the following argument. We go back to our start- 
ing point and think of the eigenangles di € [0, 27r) of 
5 G Al jv as a strongly correlated gas of N particles with 
repulsive long-range forces. (These are the particles of 
the A = 2 Calogero-Sutherland model, except that now 
they are confined to a ring of circumference 27r instead of 
L.) An important consequence of the strong correlations 
is that the total system size becomes an irrelevant length 
scale of the gas in the large- A^, or high-density, limit. 
This comes about because the (impenetrable) particles 
of the correlated gas never stray far from their average 



positions and therefore do not sense the periodic bound- 
ary conditions. There is, then, just one length scale in 
this problem for large A^, which is the mean distance. A, 
between nearest neighbors. Similarly, there is only one 
time scale: the time it takes for the wave packet of a free 
quantum particle with kinetic energy /dO^ to spread 
over a distance A. The interaction, owing to its scale- 
invariant 1 /r^ short-distance form, does not introduce an 
extra time scale. For these reasons, and because correla- 
tions always decay, the dynamical correlation functions of 
the particle gas will go to zero on space and time scales 
large compared to A and A^, respectively. Since the 
spacing of A^ particles uniformly distributed on the unit 
circle is A = 2ti/N , we conclude that the dynamical cor- 
relation functions are fully determined by the physics at 
short distances (of order 1/A^) and short times (of order 
1/N'^). Note that this is an indication of universality. 

Returning to the equivalent problem of free-particle 
motion on the symmetric space Mn, we infer that the 
curvature of A^jv becomes an irrelevant feature in the 
large-A^ limit. In concrete terms, for the purpose of 
computing the dynamical correlations near some point 
p = e'^ e AIat, we may as well replace AJa? by its tan- 
gent space Tp Al n at p. In this way the problem of free- 
particle motion on a curved space gets simplified to free- 
particle motion on a flat space. There is one complication 
however. Its nature, and remedy, are again most clearly 
seen in the particle gas picture. If the unit circle S^ is 
replaced by its tangent space, which is the real line R, 
there is nothing that holds the gas together to produce 
a finite ground state density. Fortunately, we can eas- 
ily remedy the situation by simply enclosing the gas in a 
box or any kind of confining well. More precisely speak- 
ing, what we are going to do is the following. If particles 
were added to the confining well with the size of the well 
kept fixed, the particle density would increase. There- 
fore, to have a well-defined limit when sending N ~* oo, 
we rescale the confining potential so as to maintain a con- 
stant particle density near the center of the well. Then, 
in the limit A^ — > oo a translationally invariant region will 
form around this center, the system will become ignorant 
of the confinement scale, and the dynamical correlations 
will approach the universal form we are trying to com- 
pute. 

It is now clear how we should proceed. All points of 
a symmetric space are equivalent, so without loss we 
take p to be the origin, p = =: a. We linearize 
S — expiQ around S = a, with iQ parameterizing the 
tangent space TqA^at. From the defining equations of 
Mn (i.e. = S'^ and S = CS'^C-^) we infer the 
linear conditions 



Q = Qt ^ CQ^C- 



(14) 



The correct physical dimension of length is restored by 
replacing Q {L/2n)Q. The expression for the particle- 
hole function Oxy then takes the form 
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[y-x + ieY ^y^^" 



GUQ) 



2 X — ie — Q 



Det ■ 



■ le - 



Q 



le — 



(15) 



To obtain the Green's functions of the CSM at A = 2 
and for N oo, we have to solve a one-particle problem 
with Hamiltonian 

where Cq is the Laplacian of the Euclidean space TqM n 
with metric i TrdQdQ, and ViQ) = V{kQk-^) IS a con- 
fining radial potential with a global minimum at Q = 0. 
If (Q|0) is the ground state wavefunction of H and 

K{Q,Q'-t) = {0\Q){Q\e-'''"/^m{Q'\0), 
we need to calculate the double integral 



dQ / dQ' 0^y{Q)K{Q,Q';t)0^,,{Q') 



(16) 



where dQ is a Euclidean measure for To-^n- We take 
the thermodynamic limit N oo while rescaling V in 
such a way that the density po — (0|i TtS{Q)\0) remains 
fixed. By the argument given above, the expression ( p^ ) 
with X, y, w, z at a distance of order from the mini- 
mum of the confining well, will then converge to the CSM 
particle-hole Green's function Gxywz{t) in this very limit. 

For technical convenience we choose the potential to 
be harmonic: V{Q) = mu^ Tt / Ah. We then get a 
harmonic oscillator with frequency uj, ground state wave- 
function {Q\0) = const X exp(— ma;Tr(3^/4?l) and prop- 
agator 



2Trihamujt 



X exp 



Ah sin ujt 



Tr((Q2 + Q'2) cos tot - 2QQ') 



When time is analytically continued to the imaginary 
axis, this propagator turns into the diffusion kernel for 
Brownian motion in what is called the Gaussian Sym- 
plectic Ensemble (GSE) of random matrix theory. For 
future reference we note that the basic time-ordered har- 
monic oscillator correlation function for self-dual sources 
A and B is 

{0\TTTBQ{t)TTAQ{t')\0) = ^e-i^\t-t'\ TtAB. (17) 

The average density of the GSE is known to be of semi- 
circular shape fl and in the center is given by 



po = ■\/ Nmu/hn'^ , 

so we are going to take the thermodynamic limit holding 
the product Nui fixed. 

All of the above steps are easily transcribed to the cases 
A = 1 and 1/2. For A = 1 we end up with a harmonic 
oscillator model of complex hermitean matrices Q, while 
for A = 1/2 the matrices Q are real symmetric. 



VII. SUPERSYMMETRIC CALCULATION OF 
GREEN'S FUNCTIONS FOR X = 2 



We have seen that the particle-hole Green's function 
Gxywz{t) of the CSM at couphngs A = 1/2, 1, and 2 can 
be computed as a harmonic oscillator correlation function 

Gxyn,.{t-t') = {{)\rOxy{Q{t))0^z{Q{t')m 



in the thermodynamic limit. For A = 2, the harmonic 
oscillator variable Q is a self-dual hermitean 27V x 27V 
matrix, and by using (|l^) we get 



Gxywz{t^ 



(2i£2/7r)2 



jxywz 



[t t') = (0|TC-(Q(0)C- (Q(t'))|0). (18) 



(a;-y)2(w-z)2 



Similar formulas hold for A = 1/2 and 1. Eqs. (|lq ) 
express the particle- hole Green's function as a sum of 
four terms. In Sects. VII A - VII H we calculate g'^ , from 
which g is obtained by complex conjugation followed 
by the time reversal operation. The modifications neces- 
sary for g"*""*" and g'" are described in Sect. VII 1. The 



changes that occur for A = 1/2 and 1 are sketched briefly 
in Sect. IVHI 



A. Gaussian superintegrals 

Because C^yiQ) depends on Q in a rather complicated 
way, the problem of calculating the Green's function 
is not immediately tractable. As a preparatory step, we 
are now going to make Q appear in the argument of an 
exponential function. This will reduce the problem to 
calculation of Gaussian integrals. 

Let Vka and Xka {k = 1,...,7V;q; =t,i) be a set of 
complex commuting and anticommuting (or Grassmann) 
variables. Our starting point are the formulas 

Detix±^e-Qr^^|dvd,e.^±^H^±^e^Q)v, 
Det(y ^Q) = J dxdx exp ixiv - Q)x, 

where dvdv and dxdx are suitably normalized (Eu- 
clidean) measures, and 

v{x ±is - Q)v = Vka{{x ± ie)5ki5ap - Qkaj0)vip, 

Xiy - Q)X = Xka{ySklSal3 ~ Qka,ll3)Xl[)- 

Summation over repeated indices is implied. 
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In later subsections a plethora of tensor products will 
appear, so to avoid writing a lot of indices we are now 
going to switch to basis-free notation. This we do by in- 
terpreting the integration variables Vka and Xka as the 
matrix elements of a (1-1-1) x 2Ar-dimensional super- 
matrix W. More formally, we proceed as follows. We 
introduce (canonical) bases {e'^}fc=i,...,Af, {s°'}a=^,i, and 
{e^,e^}, of particle-coordinate space C^, "spin" space 
C^, and super (or Boson- Fermion) space 14f = C^l^. El- 
ements of the dual bases are denoted by 9'' etc. We then 
take W : (^C"^ Vbf to be the linear function defined 
by 

(g> C2 by 



w 



>e Vka 



The introduction of W and its adjoint W allows us to 
combine terms: 

VkaQka,lf3Vl0 + XkaQka.l/SXiP = TtWWQ'^. 

Now let {E^'^}f^^^=B,F be an adapted basis of matri- 
ces in superspace that satisfy i?^''e'^ — e^S'^^. Defining 
{x± ie)E^^ + yE^^ we have 



W 



Similarly, we define W : Vbf 



Det- 



y-Q 



X ± ie — Q 

By differentiating both sides with respect to x and then setting y 
Tr. ' 



x±ie - 



{x ± ie)vkaVka + yXifiXW = TrVFgVF, 
which yields the formula 
= I dWdWeyi^±iT:YW{qW -WQ^). (19) 
X ± ie, we obtain 

TJ / dWdW Tr WE^^W exp ±i Tr W{{x ± ie)W - WQ^). (20) 



B. Implementing self-duality 



Formulas ( |19| ) and (^0|) allow us to write the particle- 
hole function C^y{Q{t)) as a product of Gaussian inte- 
grals over two supermatrices W^^ and W^*". Doing the 
same for its complex conjugate C~^{Q{t')), we can carry 
out the integration over Q's in the expression for the 
Green's function by setting 



Rt 



B = iW^^'w"^*' 



and using the harmonic oscillator identity 
(0| exp(Tr AQ(t)) exp(Tr BQ{t'))\0) 
= expi(0|(TrAQ(i)+TrBQ(t'))'|0). (21) 

Unfortunately, since A and B fail to be self-dual, we 
are not permitted to use the simple correlator (|l^) but 
must work with a more complicated version thereof. 
As a result, the Hubbard-Stratonovitch transformation 



(Sect. VII D) is unnecessarily complicated and the or- 
thosymplectic symmetry of the extre mal su rface emerg- 
ing in the thermodynamic limit (Sect. VII E) is obscured. 
To improve the situation wc make the following modifi- 
cation pof . 

The idea is to enlarge the supermatrix W in such a 
way that P := WW does satisfy the self-duality condi- 
tion P = CP'^C-^. To that end, we introduce an addi- 
tional degree of freedom, called "quasispin" , which takes 
two values and is thus formally identical to a spin 1/2. 
Quasispin space is denoted by Vd- (It is the cooperon- 
diffuson freedom of mesoscopic conductor physics.) W is 



extended to a supermatrix W : Vbf ® Vcd by 

adding an extra row index running through two values 
(for the two-dimensional space T4d — C^). 

Let W"^ denote the supertranspose of the supermatrix 
W. Note that W'^'^ = aW and W'^'^ = W<j, where 
a := (E^^ — E^^) (g) 1 stands for superparity. (These 
double transposition rules are implied by the relation 
(W^VF)""" = M^"'"IF"'", which is needed for reasonable su- 
percalculus.) The product WW will be self-dual if we 
impose the conditions 



W = CW^T~'^, 



W 



tw'^c-K 



(22) 



with T some c-number matrix. 



lyT ^ c-iTl^T = C-i^I^TT^T^ Since 



Consistency requires 
-C, it fol- 
lows that T must satisfy r = —ar'^. The choice we make 
is r = E^^ (g) icjy + E^'^ (g) 1. 

The extension to quasispin space in combination with 
imposition of the constraint W — CW^t^^, leaves the 
number of independent components of W unchanged. It 
is easy to see by inspection that formulas (|l^) and ( |20| ) 
remain valid if we modify W in the manner described, 
extend q to quasispin space in the trivial way {q ^ q(gl) 
and place an extra factor of 1/2 in the exponent. 



C. Integral Representation of 

We now assemble the various factors needed for the 
Gaussian integral representation of the Green's function 
g+- by (^, (^, and @). For the trace of the 
resolvent and the ratio of determinants at time t, we in- 
troduce two supermatrices W^* and W^*, respectively. 
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Similarly, at time t' we introduce a pair of supermatrices 
W^^*' and W^^*' . 

Since the multitude of supermatrices involved leads to 
lengthy expressions, we further compactify our notation. 
We do this by arranging all supermatrices vertically as a 
single giant supermatrix, denoted again by W. In tensor- 
product notation. 



Vbf (E) Vcd (E) Kr E) Vt, 



where Vt ^ ~ Kr and the subscripts "T" and "ar" 
refer to the two-Time and advanced-retarded (±ze) struc- 
tures. The new indices take the values p = t,t' and 
X = A, R. To assemble all factors, we define 



q=[{x- ie){E^^ + E^^) ® E^"^ ® -B" + 

{{x + ie)E^^ + {y + ie)E^^) ® E^'^ ® E'' ■ 



{{w - ie)E^^ + {z- ie)E^^) 
{w + ie){E''^ +E^P)i 

and we introduce the matrices 



AA 



E 



E 



t't' 



Led, 



I = E 

I' ^E 



BB 



BB 



led ' 
led ' 



AA 



' E^*, TTt = Ibfxedxar ^ E** , 



RR . 



1 



bfxedxar 



E 



t't' 



(Here and in the following we occasionally take liberty to 
change the order of the tensor products, for notational 
simplicity. What is meant will always be clear from the 
index nomenclature used.) With these definitions, 
has the expression 



9xywzi^ < ) — 2 



dWdW X 



Tt{WIW) Tr(W^/'M^)e 
(0|r exp § Tr VFr/(7rtM^g(t) 



iTrWriqW/2 



Trt'WQ{t'))\Q). (23) 



A diagonal matrix rj containing elements ±1 was inserted 
for convergence of the bosonic integrations, see ( |l9|) and 
(pO|). The sign of these matrix elements is fixed on the 
bosonic subspace but is arbitrary on the fermionic one. 
The "good" choice for rj turns out to be [ppl|] 



BB 



E 



FF , 



lar) ® ledxT- 



We make this choice because it leads without further ado 



to the metric of the extremal surface of Sect. |VII E| being 
Riemannian (as opposed to being indefinite). In contrast, 
if we made the more natural choice rj' = a^^ (E IbfxcdxT, 
we would eventually be forced to rotate integration con- 
tours to achieve positive definiteness. 

We add a brief discussion of symmetries. Recall that 
W and W are related by W = CW^t'^ with 



+ i;^^®lcd)(»larxT. 



For t = t' and q a multiple of the unit matrix, the expo- 
nential part of the integrand in (p3|) depends on W only 



through the combination WrjW, which is invariant un- 
der transformations W gW, W Wg'' if g satisfies 
7]. Invariance of the constraint W = CW^t~^ 
The two conditions on g can be 
combined into the equations 



requires g^ — Tg"^T~^ 



g = {Tri)g 1 (tijY 



-it - 
rjg rj 



(24) 



The first of these defines an orthosymplectic complex Lie 
supergroup Osp(8|8). The second fixes a noncompact 
(pseudo-unitary) real subgroup Uosp(4, 4|8), which we 
shall denote by G for short. G is the symmetry group of 
our problem and will play a central role later. 



D. Hubbard-Stratonovitch transformation 

Formula (|2^ ) is useful because it enables us to inte- 
grate out the matrix variables Q and W and be left with 
a small number (independent of N) of integrals to do. 
This is achieved as follows. 

Combination of the harmonic oscillator identities (|Tl 
and (|2^) gives 

(0|rexp(Tryig(t))exp(TrBQ(i'))|0) 
= exp ?iTr(yl2 + + 2ABe~'"l*"*'l)/2TOtj. (25) 

We apply this relation to (^3|) by identifying 

iWrjT:tW/2^ A, iWrjT:t'W/2 = B. 

Putting a := i'KtWWrj/2 and b := iTTt'WWrj/2 and using 
the cyclic invariance of the trace, we rewrite the right- 
hand side of (E^) as 



exp h STr(a2 + 2a6e~''^l*^*'l ) /2mw. 



(26) 



The symbol STr denotes the supertrace. The next step 
is to linearize the quadratic form + fo^ + ... by the in- 
troduction of a Hubbard-Stratonovitch field M . Let M 
be a supermatrix which has the same symmetries as the 
composite object WWrj, and let a Gaussian Af-average 
be defined by 

= J dM • exp-mujV{M)/2h, 

V{M) = STr(M2 + + 2e-l*-*'lM,,, Af,, J, (27) 

where Mij = niMnj for i^j = t,t'. We then claim that 
the expression (M) is equal to 



(exp STr M{a + b)) m ^ (exp f Tr WrjMW) 



M- 



Formal verification of this claim is the simple matter of 
doing a Gaussian integral by completing the square. On 
a rigorous level, however, one needs to demonstrate that 
the integration contours for the bosonic-bosonic variables 
in M can be chosen to be uniformly convergent in W. 
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This is a tricky matter which will not be discussed here 
but is well understood; see Ref. H^, Sect. 5.2. 

Aside from being Gaussian in W, the last expression 
has the crucial property of being invariant under unitary 
transformations W WU, W WW [U e Sp(27V) C 
U(2A^)). This means that the particle-coordinate and 
spin degrees of freedom have now been completely de- 
coupled. Hence, doing the Gaussian integral over the 
supermatrix W subject to the constraint W = CW^t~^, 
we get the inverse of a superdeterminant raised to 
power: 



the iV'h 



dWdW e'TrW^iM-g)w/2 ^ sDet-^(M - q). 



By differentiating this equation with respect to the pa- 
rameters {x — is) and {w + ie) which are contained in 
q and couple to Ti W rj I W and TrWrjI'W, we can re- 



express g"*" in (ESp as an Af-average: 



Z{M) = STr 



Z{M))m 
I 



STr- 



SDet-^(M-9). (28) 



M -q M -q 



Note the close similarity to the starting expression ( |lq ) 
and ( p^ for : the resolvents and determinants in 
Q-space have simply been replaced by resolvents and su- 
perdeterminants in Af-space. 



E. Thermodynamic limit 

The Gaussian integral transformation from Q to M has 
reduced the number of integration variables enormously, 
from order N"^ to order . The large number iV now 
appears explicitly in the integrand, allowing us to make a 
saddle-point approximation which becomes exact in the 
thermodynamic limit N oo. 

Recall th at we take N oo holding the particle den- 
sity po = y/WnuJjTm^ fixed. We anticipate that the su- 
permatrix M will be of order at and near the saddle 
point, whereas q (containing the positions x, y, w, z) and 
|t — i'l are of order N^. To extract the leading large- iV 
behavior of the integral ( p8| ) and (|2^) , we look for those 
M where the superfunction 

n{M) = SDet~^(Af)exp-mwSTrAi"2/2?i. (29) 

is maximal. By the relation 

SDet"^(Af) = exp-TVSTrlnAf, 

such supermatrices are solutions of the saddle-point equa- 
tion NM"^ + mLoM/h — or, equivalently, 



Af" 



iN/k^) 



with fco = Trpo- They form supermanifolds that are in- 
variant under "rotations" M ^ gMg~^ by elements g of 
the symmetry group G. 



Let us now make M dimensionless by substituting for 
it the rescaled expression {iN/ko)M. The solution spaces 
of the rescaled saddle-point equation AP — 1 can be la- 
belled by the eigenvalues of M, which are equal to ei- 
ther -1-1 or —1. It turns out that such eigenvalues do 
not lie on the integration contour that is suggested by 
the requirement of uniform convergence for the bosonic- 
bosonic variables. However, a subset of the surfaces on 
which ri(Af) is maximal can be reached by deformation 
of the integration contour using Cauchy's theorem. As 
far as the bosonic-bosonic sector is concerned, the signs 
of the eigenvalues of M on such a "reachable" surface 
are uniquely determined by the pole structure of n{M), 
which in turn is determined by the small imaginary parts 
±ie, omitted in (p9|); see |2|], Sect. 5.3. As for the 
eigenvalues of M associated with the fermionic-fermionic 
sector, no constraints from convergence and analyticity 
exist. In this case, the determining agent is the limit 
iV — > oo, which singles out the extremal surface with the 
largest dimension d — db ~ d{ (bosonic dimension minus 
fermionic dimension). By inspection one finds that d is 
largest for the surface containing the element 



A = af 



1 



bfxcdxT- 



We restrict the integration in ( |28| ) to this surface by set- 
ting M = gAg-^. Expanding SDet"^(i7Vfc(7^Af - q) 
with respect to q and exjpiLu\t — t'\ with respect to |t — 
and keeping only the terms that survive in the thermo- 
dynamic limit, we find 

g+- = -(fco/2)2 / dfi{M)STT{IM)STr{I'M) 



X exp i^iko STr qM + iuQ\t - t'\ STr Mn'Mft) (30) 

with ujQ — hk^lm. As usual, integration over the Gaus- 
sian fluctuations around the surface produces the uniform 
(or G-invariant) integration measure dfi{M). 



F. Choice of polar coordinates 

What we did up to now was standard technology 
and obvious to the expert. Now the real work 
begins! Our task is to calculate the integral (|30|). We 
begin with some notational and conceptual preparations. 

The supermatrix M has dimension 16 x 16, since the 
spaces Vbf, Vcd, Kir and Vr, on the tensor product of 
which M acts, are all two-dimensional. Recall the defi- 
nition of the symmetry group G as the set of 16 x 16 su- 
permatrices g that obey the conditions (|2j). We endow 
G with its natural metric structure given by STr dgdg~^. 
Let K be the subgroup of elements fc € G that commute 
with A. Since M = gAg~^ is invariant under g <—>■ gk 
(fc e K), the integrand in ( ^0| ) is a function on the coset 
space G/K. Elements of this coset space are denoted by 
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gK, and the uniform integration measure on G/K is de- 
noted by dgK- Thus, the integral (|30| ) is of the general 
form 

/ f{g^)d9K. (31) 
Jg/k 

To calculate such an integral we need to choose a suitable 
system of coordinates, or parameterization, of G/K. In 
making this choice we are guided by the general form 
the particle-hole Green's function is expected to have: a 
good coordinate system should contain the velocities of 
the elementary excitations produced by the particle-hole 
operator as a subset, and it should make the integral 
separate into two form factors corresponding to the final 
time t and initial time t' . In the present subsection and 
the one that follows, we will describe the parameteriza- 
tion that meets these expectations. 
Let T3 be defined by 

T3 = TTf - TTf = Ibfxcdxar ® 0"^. 

We denote by Go the subgroup of G whose elements 
commute with T3. In other words, Gc is that part of 
the symmetry group which operates at fixed time. Fur- 
thermore, we denote by A the six-dimensional abelian 
subgroup of G generated by (i = 1, 2) 

i/l led® (C+^Kxar, 

Hl+2i — ® El^^ (g) (C^j^)Txar, 

H2+2r = E^^ ® E';^ ® {Cj^,A)Tx.r, (32) 

with Cfj = (g) E'l ±E*'*(g) E^\ These generators are 
off- diagonal with respect to the chosen bases of Vt and 
Var, which implies that they anticommute with and A. 
Therefore, the elements a of the group A satisfy 

Technically speaking, A is a maximal abelian subgroup 
for the Cartan decomposition of G w.r.t. K. 

Given Ge and A, we introduce nonstandard polar co- 
ordinates on G/K by the map 

(j) : Ge/M X A+^ G/K, 
igM,a)^ gaK, 

where M is the subgroup of elements of Kc = K n Gc 
which commute with all elements of A, and A+ is some 
connected open subset of A such that cj) is bijective. The 
abelian factor a contains the degrees of freedom which 
are "radial" with respect to the polar coordinate decom- 
position, while g € Gc will be parameterized by a set of 
"angular coordinates" . 

By the substitution rule for superintegrals , the po- 
lar coordinate map transforms the invariant integral 
(|l|) into 



{L ^^^"^^'^^) + - ^^^^ 

Here dg is a (suitably normalized) Haar-Berezin measure 
for Gc, da is a Euclidean measure on A, and J is the 
superjacobian of the transformation. The dots indicate 
correction terms (so-called "boundary terms"), which are 
due to the super nature of the integral and will be dis- 
cussed later. 

The jacobian J is calculated in Appendix A. If we put 
Ina = X^Li ^i^i + 'PiHi+2 and introduce 

vi = —Ws cosh 201, V2 — cosh 26*2, 

Wl = — cos 2(y9i, V3 = Vs cos2(p2, 

V2 — —Vg COs2(p3, V4 = WsCOs2(y94, (34) 

our result for J = li^C)^ normalization) is given 

by (H) with A = 2, rip = 2 and rih = 4. We choose for 
A+ the fundamental positive domain defined by 

— OO < Vi < —Vs < Vi < V2 < V3 < V4 < Vs < V2 < OO. 

We will now extract from the argument of the exponen- 
tial function in ( pO| ) those parts which depend solely on 
the abelian factor a parameterized by {vi, Vj}. The term 
multiplying |t — t'| is evaluated as follows: 

^iSTrMu'Mt't = - STr(l + r3)M(l - T3)M 
= STrr3MT3Af = STr (r3aAa"^)^ = STra^ 
= STr cosh 4 In a = 2 STr cosh^ 2 In a 

= 4 E - 4 E ^1 =■■ ^(^Eia)/mvl (35) 
^ i ^ j 

The term STr qM containing the coordinates x, y, w, z of 
the points of particle creation and annihilation is rewrit- 
ten as 

STrqM ^ STr q{ga^ Ag-^) 
= STrqiga'^Ag-^ - a^A) + STr Ag cosh 2 In a. 

Defining mh = —m/2 and introducing 

Pp = mivi + V2) + TOh(wi -I- V2), Ph = m^{vz H- W4), 

we can write the ^-independent term in the form 

—mVa STr A(7 cosh 2 In a = (x — w)Pp(a) + (y — z)P\^{a). 



G. Identification of the form factor 

The introduction of polar coordinates on G/K by the 
map (j) separates the integral (^^ over M into a (radial) 
integral over the abelian group A and an (angular) in- 
tegral over the symmetry group Go operating at fixed 
time. The latter, in turn, separates into two commuting 
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factors, Go = x Gf, one for each time. Defining for tical, let us concentrate on the first one. We introduce 



i = t,t' the projections 

Qi — TTiQ, Ki{a) = fco7ri;Acosh21na, 

we get 



At = TTtA and ttk = E^^ ®E^^'®lcd®E^\ TTp = Trj-TTh. 
Then, 

qt = xTTp + j/TTh - ieAt 

= ^{x + y)TTt + 5(2; - v){t^p - TTh) - ieAt, 



Hence, the integral over Ge splits into two factors, one and, since STry ^i^tg — STrTTt = for g £ Gt, the inte- 
for Gf and another for Gf ' . Since these are formally iden- gral over Gf is given by 



fl{x~y-a)^ / STr(/5A'f(a)g-i) exp - iSTv{\{x - y){T:,, - -Ky,) - ieKt){gKt{a)g-^ ~ Kt{a)) 

JGt '- 



(36) 



To summarize, our complete expression for g^ is 

atywzit)^ / /f^(a;-2/;a)/^(w-z;a)exp[-it£:(a)/n + i(x-u;)Pp(a)/;i + «(j/-z)Ph(a)/?i]J(a)da. (37) 



Having gone through a long and indirect derivation, 
we are finally in a position to interpret the mathematical 
structures at hand and translate them into physics. We 
see that a particle-hole excitation of the ground state of 
the Calogero-Sutherland model at A = 2 fractionalizes 
into six elementary excitations. The velocities of these 
excitations appear as the degrees of freedom of a six- 
dimensional abelian group A in our treatment. There 
are two particle excitations with velocities vi , V2 and 
four hole excitations with velocities ui, ...,-1)4. Note that 
Vi < —Vs and V2 > +Vs are on opposite sides of the 
pseudo-Fermi sea. The radial functions 

P{a) — hmVs STr7rf/Acosh21na, 



E{a) 



'STrcosh2 21na, 



are the total momentum and energy of the particle-hole 
excitation. The equation P = Pp + Py, decomposes the 
momentum into its particle and hole components. Our 
expressions for Ph and Pp show that a single hole frac- 
tionalizes into two hole excitations (with velocities W3, U4) 
whereas a single particle fractionalizes into two particle 
and two hole excitations (with velocities wi, W2; wi, ^2). 
Finally, 



F{x;< 



= — 2ftix; a)y/ J {a) 



(38) 



is identified as the form factor: 



— ^(2; _ y- Q-)gi3:Pp(a)/R+iyPh(a)/fi_ 

The last statement is not affected by the fact that g^ is 
only one out of four contributions <?'^^ (ar = ±1) to the 
particle-hole Green's function. The other contributions 



yield the same form factor; see Sect. VHI. 

It is possible, though by no means easy, to calculate 
the form factor F in closed analytic form by doing the 



integral in (|3^). The result is neither simple nor illumi- 
nating, and we will not give it here. Substantial sim- 
plifications occur for |a; — ?/| ^ 1/fco- In this limit, 
F{x — y; a) ~ Pp(a) x Ph(a) becomes independent oix — y 
and separates into two factors for a single particle and a 
single hole. The technical steps of taking |x — ?/| 00, 
e ^ 0, and carrying out the final integrations, are done 
in Appendix B. Our result for the single-particle form 
factor has been presented in Sect. IV. The result for the 



single-hole form factor is the one given previously 



H. Boundary terms 

Berezin's proof p2| of the substitution rule for superin- 
tegrals goes through without modification only for func- 
tions with compact support, or if the integration domain 
has no boundary. For polar-coordinate superintegrals 
such as (^3|) , whose radial domain A+ does have a bound- 
ary, one expects extra terms (indicated by dots) to ap- 
pear in general. The structure of such "boundary terms" 
was first investigated systematically by the mathemati- 
cian Rothstein |^^. Unfortunately, Rothstein's general 
theory is too implicit to be useful for anything but the 
simplest applications. In the present problem we are 
dealing with the polar-coordinate integral on a super- 
manifold with the rather special structure of a symmetric 
space. For this powerful method for construction 

of all boundary terms has recently been developed [^ . 
The basic idea underlying this method is sketched for a 
simple example in Appendix C. Postponing the full ex- 
position of the method to a future publication, we are 
now going to state the result obtained. 

It turns out that all the boundary terms that occur for 
polar-coordinate superintegrals of the type (|3^), are as- 
sociated with nonintegrable singularities of the radial in- 
tegration measure J{a)da. While such singularities never 
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exist for ordinary (i.e. nonsuper) polar-coordinate inte- 
grals, they do exist in the present case because the su- 
perjacobian, being a superdeterminant, puts factors both 
in the numerator and in the denominator. By the na- 
ture of the polar coordinate map, the singularities all lie 
on the boundary 9A+ of the radial space A+. For our 
choice of A+, nonintegrable singularities occur for (1) 
vi = vi = V2 = -Vs, and (2) U3 = U4 = U2 = +Vs. 
Each of these conditions defines a three-dimensional sub- 
space {i = 1,2) of A+. Their intersection defines 
a zero-dimensional subspace, which is the unit element 
a — 1. Correspondingly, there are three boundary terms 
in the present problem. To formulate them, let be 
the subgroup of elements of Kc which commute with all 
elements of A^. Furthermore, let Ji : A^ R be the 
reduced jacobian which is obtained by omitting from J 
all singular factors. With these definitions, the complete 
formula for transformation to polar coordinates is 



figK)dgK = / figaK)dgM J{a)da 

G/K J A+ yjGe/M J 

+ J2 f . i I /(5aK)dgM. I Ha)da 



i=l \JGe/Mi 



Go/Ke 



fig'^)dgK^- 



When the function / is given by the integrand of (30), 
this formula is interpreted as follows. The first term on 
the right-hand side (the "regular" term) corresponds to 
fractionalization of V'lV'yiO) into six elementary excita- 
tions, as discussed above. What the appearance of two 
boundary integrals with diniA^ = 3 means is that there 
exists the additional possibility for V'lV'ylO) to fraction- 
alize into three elementary excitations. Two of these are 
holes, and one is a particle. (It is not hard to show that 
the regular term vanishes by symmetry for x — y, so that 
the dynamical density correlation function is completely 
determined by fractionalization into three elementary ex- 
citations, in agreement with the findings of Q.) Finally, 
the last boundary term, where the abelian factor a is set 
to unity - i.e. no elementary excitations are present - 
accounts for the fact that the ground state is contained 
in ipl,ipy\0) with nonvanishing amplitude in general. In 
other words, this term is equal to the product of ground 
state expectation values (O|'!/'J,'0a;|O)(O|7/'t^Vz|O); which de- 
termines the long-time limit of the dynamical correlation 
function. 

All of the boundary terms make nonvanishing con- 
tributions to g'^ , and hence to the full particle- hole 
Green's function, in general. However, when the points 
X and y (or w and z) are removed from each other, only 
one of the boundary terms survives. This is the one that 
results from setting vi — vi — V2 = —Vs- To extract its 
contribution to the advanced part of the single-particle 
Green's function, we proceed as in Appendix B. The re- 
duction in dimensionality leads to somewhat simplified 



calculations in comparison with the regular term. The 
form factor is now a function of just a single variable 
V ^ V2. Doing the calculation, and including the con- 
tribution from g ^, we obtain the result for the elemen- 
tary term in the single-particle Green's function given in 



Sect. IV 



I. The terms g^"*" and g 
Formula ( p^ ) for the particlc-holc Green's function 

involves four terms: g^^, g^ , g ^, and g . We have 
shown in detail how to calculate g^ , which is related to 
g by complex conjugation followed by time reversal. 
The calculation of g is identical except for the minor 
changes we are now going to indicate. 

g differs from g"' only by C^y{Q{t)) replacing 
Cly{Q{t)) in (|l|). Consequently, to get g— from (H), 
all we need to do is to replace the definitions of / and q 

by 



/ = E""" (g led ® E" (g) E", 

q = [{{x - ie)E^^ + {y- ie)E^^) (g, E^^ (g, E'' -f 
{x + ie){E''^ + E^^) g) £;«^^ g) + ...] g) led- 

The supertrace in the expression for the phase factor 
exp —iko STr qA cosh 2 In a changes to 

—Vg STr gA cosh 2 In a — 2{x — w){vi + V2) 
-{y - w){vi + V2) - {x- z){v3 + V4). 

This equation illuminates the physical content of g : 
since the variables Vi and vj are the velocities of the el- 
ementary excitations present in iplj{t')^pz{t')\0) , we see 
that we are now dealing with a situation where the two 
holes created by ipz{t') are annihilated by tpx{t) (rather 
than by tpl{t) as before), while the two holes created by 
tpi,{t') are annihilated by ipl{t). It is clear from this pic- 
ture, and it is not hard to verify by direct calculation, 
that g vanishes in the limit |a: — ?/| 00, x — w and 
y — z held fixed. The same statement applies to g^^ ■ 
Therefore, the single-particle Green's function receives 
no contribution from g and g"*"*" but is solely deter- 
mined by g^ and g ^ . 



J. The cases A = 1/2 and 1 

A major benefit from our basis-free approach is the 
ease with which the transcription to the case A = 1/2 can 
be made. We now briefly describe the essential modifi- 
cations, paying no attention to changes of normalization 
factors etc. 
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The harmonic oscillator variable Q becomes a real sym- 
metric matrix of dimension N x N, while the basic time- 
ordered harmonic oscillator correlation function ( p7| ) re- 
mains formally the same (for symmetric A and B). The 
expression for the particle-hole function Oxy is replaced 
by 



CUQ) = Tr 



y/y~ X + ie 
X ~ ie — Q 



'Det 



y + le ■ 



X + ie — Q 



The trace of the resolvent is treated in the same way 
as before. To implement the symmetry requirement 
(W^VF)"^ = WW, we take to be a hnear function 
W : — > Vbf (E) Vcd subject to the constraint 

W^tW^, W = W'^t-\ 

Self-consistency now implies r = +(7t'^, and we choose 



E 



BB 



l + E 



FF 



Note that this is different from the case A = 2 only by the 
exchange of the bosonic-bosonic and fermionic-fermionic 
sectors. Consequently, the cases A = 2 and A = 1/2 
are connected by a duality transformation that simply 
exchanges these sectors. (See (p^).) 

To deal with the square root of the determinant, we 
write it as 

Det(7; + ie~Q) 



^J'Det{y + ie-Q) Det{x + ie - Q) 

For the numerator we introduce a complex fermion, and 
for the denominator two real bosons. (This is the same 
count as in the advanced (— ie) sector.) 

The Green's function defined in ( p^ ) is still given 
by (HI) but with 



q^[{x - ie)lbfxcd ' 
{y + ie)iE^^(> 



{x + ie)E 



BB , 



EE 

eI\ 



AA 



E 



FF , 



E 



RR-i 



' led) «) E 



RR 



The dots indicate analogous terms corresponding to time 
t', which are unchanged. Hubbard-Stratonovitch trans- 
formation and thermodynamic limit remain the same. 

The symmetry group G is defined by ( p4[ ) (with the 
modified expression for t), and the abelian group A ap- 
pearing in the polar decomposition of G is again six- 
dimensional. From the duality symmetry connecting 
A = 1/2 with A = 2, we readily see that four of the 
six degrees of freedom of A are now particles and two 
are holes. We choose the fundamental positive domain 
A+ 

— OO < Vi < V2 < —Vs < Vi < V2 < Vs < V3 < V4 < OO. 



The decomposition of total momentum P into the com- 
ponents Ph (coupling to y — z) and Pp (coupling to x~w) 
now reads 

Ph = m{vi - 2^2), Pp = ■m{vi +V2 + V3- 2vi). 

From this we learn that a single particle removed from 
the ground state fractionalizes into a particle (with mass 
to) and a hole excitation (with mass — 2to), whereas a 
single particle added to the ground state fractionalizes 
into three particles and one hole excitation. 

The expression for the particle-hole form factor is re- 
placed by 



F{x;a)^ lim J-Jt{x;a)yJ J{a) 



where /f is still given by (36) but now 



{E 



BB 



E 



FF 



1 



cdj 



E 



RR 



by the changed formula for q. Evaluation of the single- 
particle form factor proceeds in essentially the same way 
as in Appendix B. 

Finally, the case A = 1 (free fermions) is trivial in the 
CSM but nontrivial in the corresponding random matrix 
or harmonic oscillator formulation. Thus the balance is 
now reversed: rather than solving the former using the 
latter, we can put this case to use by building up our 
superanalytic skills from the free-fermion results. The 
symmetry group of the supersymmetric integral repre- 
sentation of the particle-hole Green's function at A = 1 
is a pseudo-unitary group G = U(2,2|4) C Gl(4|4). 
The dimension of a maximal abelian subgroup A for 
the symmetric space G/K is four (two "particles" and 
two "holes"). With some effort, one can show that the 
regular term in the corresponding supersymmetric polar- 
coordinate integral vanishes. Thus the Green's function 
is fully given by the boundary terms (with two radial de- 
grees of freedom, one particle and one hole) in this case. 
This is as expected, of course, since ip^ipy acting on a 
noninteracting Fermi sea produces exactly one particle 
and one hole. 



VIII. CONCLUSION 

In conclusion, we have completed our earlier study of 
dynamical correlations of the Calogero-Suthcrland model 
by obtaining the advanced part of the CSM single- 
particle Green's function at the special couplings A = 1/2 
and 2. Our approach is based on the relation of the CSM 
to radial motion in a harmonic oscillator model, where 
the oscillator variable is a real symmetric N x N ma- 
trix (self-dual hermitean 2N x 2N matrix) for A = 1/2 
(A — 2). The relation becomes exact in the thermody- 
namic limit of infinite particle number N. 

To obtain the dynamical correlations of the harmonic 
oscillator model, we used superanalytic methods that 
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were developed in the context of random matrix theory 
and mesoscopic conductor physics. The CSM particle- 
hole Green's function for A — 1/2 and A = 2 was shown 
to be a sum of four structurally identical terms, each of 
which is an integral over a Riemannian symmetric super- 
space G/K. The spaces G/K for both couplings are real 
sections of the complex(ified) supermanifold Gc/Kc — 
Osp(8|8)/Osp(4|4) X Osp(4|4). They transform into each 
other by a kind of duality (or particle-hole) transforma- 
tion which exchanges the bosonic-bosonic and fermionic- 
fermionic sectors. The space G/K for A = 1 is a real 
section of G1(4|4)/G1(2|2) x Gl(2|2). 

The key to further progress was the identification of 
a suitable polar decomposition of G/K into a maximal 
totally geodesic flat submanifold ("radial" space) and a 
transverse manifold ("angular" space). Making this de- 
composition corresponds to expanding the particle-hole 
Green's function in a basis of energy eigenstates: 

The radial coordinates of G/K were interpreted as the 
velocities of the elementary particle and hole excitations 
which tpl,tpy\0) decomposes into when — s- oo. Thus, 
doing the radial integrals amounts to summing over i' in 
that limit. Moreover, the transverse manifold was seen to 
be a direct product of two supermanifolds, one associated 
with initial time t' and the other with final time t. Inte- 
gration over (either) one of these yields the particle-hole 
form factor (j/|V'|i/'y|0), see (H |8|). 

The radial space of G/K has both compact and non- 
compact degrees of freedom. On making the identifica- 
tion with elementary excitations of the CSM, the for- 
mer translate into holes, and the latter into particles. 
For A = 1/2 (A = 2) we find two (four) holes and four 
(two) particles. The generalization of this to rational 
X = p/q is that there will be 2p holes and 2q parti- 
cles in ViV'yiO) in general. However, this not the whole 
story yet. When the substitution rule is applied to a 
superintegral /q/^: fiO^C^ it is found that such an inte- 
gral expressed in polar coordinates, must be corrected by 
the addition of several lower-dimensional integrals (the 
"boundary" terms) . These corrections are a novel feature 
of superintegration, and are needed to restore invariance 
of the integral under translations £, g ■ (, {g & G), 
which is lost when the naive substitution rule involving 
only the "regular" term is assumed. In our case there are 
three distinct boundary terms. For A — 1/2 (A — 2), two 
of them have a radial integration domain which is three- 
dimensional, with two (one) of the radial degrees of free- 
dom being noncompact, i.e. particle-like, and one (two) 
compact or hole-like. The remaining boundary term ac- 
counts for the ground state component in V'lV'ylO)- For 
X = y the regular term in the polar-coordinate integral 
vanishes, and we recover the result that the dynamical 



density correlations at X = p/q are exhausted by contri- 
butions from p holes and q particles. The single-particle 
form factors were extracted from {vltpltjjylO) by separa- 
tion of points {\x — y\ oo). 

The methods described in this paper readily extend to 
many-parti cle G reen's functions. By following the steps 
of Sects. |v| |VIl| wc can, in principle, derive an exact ex- 
pression for the form factor of any excited state {v) which 
is produced by several actions of and on the ground 
state. Such a form factor always separates into two parts. 
The first one, F^'^\ is universal and is given by (^. The 
second one, , is given in the form of an integral rep- 
resentation analogous to (^,^8|). For example, for A 2 
we get 

F^'Hx^...x^W) = r[{x,-x,r [ e^T.^..^,-^^^ 

J0{2n) 

V = diag(wi,?;i, ...,?;„, u^), X = diag(a::i, xi, ...,a;„,a;„) 

where dg is a Haar measure for 0(2ri). (Actually, be- 
cause of the twofold degeneracies of X, it is sufficient to 
integrate over 0(2n)/0(2)"-.) 

As a final remark, we are intrigued by the finding that 
the p holes and q particles excited by the CSM density 
operator, can be interpreted as radial degrees of freedom 
of a symmetric supermanifold G/K for A = 1/2, 1, and 
2. We wonder whether there could exist a "quantum" de- 
formation (with deformation parameter A = p/q) of the 
symmetric supermanifold, such that the maximal "torus" 
of the deformed manifold had p compact and q noncom- 
pact directions. If such a deformation existed, the treat- 
ment of our paper might carry over to more general values 
of A, and the problem of calculating correlation functions 
for the Laughlin wavefunction in one dimension and in 
the thermodynamic limit, might be solved. 
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APPENDIX A: CALCULATION OF A JACOBIAN 

We are going to calculate the jacobian J associated 
with the map (/) : Gc/M x A+ ^ G/K, (.gM, a) ^ gaK. 
According to its definition by the substitution rule for 
superintegrals [Q, J is the superdeterminant of the su- 
permatrix that expresses the differential of (j) in an or- 
thonormal basis of the tangent spaces of Gc/M x A+ and 
G/K. Let Q, JC, A, Gc, and M be the Lie algebras of G, 
K, A, Gc and M, respectively, and let P be the tangent 
space of G/K at the origin 1 • K. Since the invariant 
measure on G/K is not changed by left translation, J (a) 
equals the superdeterminant of the linear transformation 
Ta:{ge-M)xA^r given by 
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as 



s=0 



where Ad{a)Z — aZa^^ denotes the adjoint action, and 
Z-p means the T'-component oi Z E Q with respect to the 
orthogonal decomposition Q = JC + V. 

To calculate J(a) — SDet Ta, we further decompose JC 
and V into their even parts, /Ce and Pe, and odd parts, 
/Co and Vo^ with respect to the involutory automorphism 
Z T3ZT3. Since the elements of ^ C T'o anticommute 
with both A and T3 , the following commutation relations 
are obvious: 



From these we see that, ii Z = X + Y is the decomposi- 
tion oi Z hy Gc- M ^Vc + (/Cc - M), then 

Ta{X + Y,H) = + coshad(lna) X - sinhad(lna) Y 

where we have introduced ad(lna)X := [lna,X] (com- 
mutator) and used the relation Ad(a) — expad(lna). It 
follows that 



J (a) — SDet coshad(lna) 



SDet ( sinhad(lna) 



The superdeterminant of a supermatrix equals the 
product of its bosonic eigenvalues divided by the product 
of its fermionic ones. We thus need the eigenvalues of the 
linear operators 

ad(ln a) : JCo + ^ + ICo, 
ad(ln a) : JC^ + Vo ^ Vo + IC^- 

Because ad(ln a) is block off-diagonal, these occur in pairs 
with opposite sign. Let us denote a set of positive eigen- 
values by and A^, respectively. Furthermore, let 
\ma\ denote the multiplicity of eigenvalue a, and put 
sign(mQ.) = -f 1 (—1) for a bosonic (fermionic). Then the 
expression for J becomes 



J(a) 



^sinhQ;(lna)^ JJ^ |^cosh/3(ln a) 



aeAJ 



Note that the derivation of this formula has used no 
more than the generic geometric structures underlying 
any symmetric space G/K with an involution T3. It is 
therefore valid for all of the cases A = l/2, A = l and 
A = 2. 

Calculation of eigenvalues a (also called roots) and 
multiplicities TOq is a standard exercise in linear alge- 
bra. Adopting the parametrization In a = X]i=i ^i^i + 
'PiHi+2 with Hi given in ( |3^ ) we find the following 
eigenvalues (resp. multiplicities jfor A = 2: 



A+ : 9i±92 (+4), 

iipk ± i(pi (+1) {k < l; k + I odd), 

9k ±m (-2) {k + l odd), 
A+ : 29k (+3) (fc = l,2), 

iipk ± j<y3; (+1) {k < I] k + l even), 

9k ± i(pi (-2) (fc + I even). 

The roots a for A = 1/2 are obtained from this by a 
simple duality transformation a ia. 

The jacobian J is easily converted to the form given in 
(^), J cx jF'^'^^p, by using the trigonometric identities 

2 smh{A + B) sinh(A - B) = cosh(2A) - cosh(2B), 
2 cosh(74 + B) cos\v{A - B) = cosh(2A) + cosh(2B), 

and recalling the substitutions (|^ . The duality symme- 
try (|l^) relating A = 2 with A = 1/2 is a consequence of 
the duality of the corresponding root systems (a <-*■ ia). 



APPENDIX B: CALCULATION OF THE 
SINGLE-PARTICLE FORM FACTOR 

We are going to extract the single-particle form factor 
_Fp(a) = {Q\tpx\o) exjp ^ixP-p{a) /h from the particle-hole 
form factor F{x — y; a) by taking jx — y| — > 00 in (|3^ ) 
and (|3^). In this limit, the integrals over some of the 
degrees of freedom of Gt — Uosp(2,2|4) can be done by 
stationary phase. We begin by isolating these degrees of 
freedom. 

First of all, we make the transformation g i— > g^^ which 
turns the term coupling to x — y into STr iirt(a)(7(7rp — 
TTh)g^^ and, being an isometry, leaves dg unchanged. We 
abbreviate TTp — TTh =: D. Next, we decompose Gt by 
Gt ~ (Gt/H) X H where 

11 = {he Gt\hDh-^ = D}, 

and we write g £ Gt as g = sh, accordingly. The fac- 
tor /i e H drops out of the expression gDg^^ = sDs~^. 
Moreover, the measures match under factorization, i.e. 
dg turns into dsndh, with dsn {dh) being the invariant 
measure on Gf/H (H). Now, because the phase of 

A{s) = exp-^{x - y)STT Ktia){sDs-^ - D) 

oscillates rapidly for |a; — y| ^ l/^Oj the integral over 
Gt/H can be done by using the stationary phase ap- 
proximation. Setting s = expAT and doing the Gaussian 
integral that results upon expansion to lowest nonvanish- 
ing order in AT, we get 



(x-y) M A{s)dsH 

JGt/H 

^ {V2 - Vi){vi - Vi) 



\x — y\ — >oo 



n 
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This expression precisely cancels that part of the phase 
space factor -^J which connects the velocities of the 
"particle set" (wi, vi, ^^2) with those of the "hole 
set" {v3,V4^), so that \/j gets replaced by the product 



4°^ X (Eq. (|)). By the particle-hole 
structure of D = TTp — TTh , the group H decomposes as 

H ~ Hp X Hh = Uosp(2, 2|2) x 0(2). 



Hence, the form factor separates as expected: F 
Fp X Fh , and we have Fp = F^"^ x F^^\ Fh = F^ ^ x Fj^^^ . 
The factor pertaining to the single hole is trivial, 



FjI^^ (a) = / dh = const, 
while the single-particle form factor Fp^'' (a) equals 

"'Hp 

where a normalization constant has been omitted. 

The remaining task is to carry out the integration over 
the group Hp. Simplifying the notation by writing K for 
Kt{a) and omitting a term eSTr AjAT, which vanishes in 
the limit e — > 0, we encounter the integral 

dh STr{KhIh-^) exp -e STrKhAth-^. 



Hp 

We will do this integral in several steps. The first 
one is to factor Hp by Hp ~ (Hp /Ha) x Ha with 
Ha — Uosp(2|2) x Sp(2) being the subgroup of Hp that 
fixes At w.r.t. the action 1— > hkth^^ . Again, the 
invariant measures match. Ha is a maximal compact 
subgroup of Hp, so Hp/HA is a (super-)symmetric space 
of the noncompact type. The integrand is constant on 
Sp(2) and depends on Uosp(2|2) only through the com- 
bination hlh~^ in the pre-exponential factor. Therefore, 
the integral over Ha is trivial and can be done by "or- 
thogonality" . (By orthogonality we mean the relation 

/ dU V{U)o,pV{U-%s = co(-l)l''l(5„5(5^^, 

JG 

valid for any irreducible representation I? of a compact 
supergroup G.) Thus, with 7r+ being the Uosp(2|2) scalar 
7r+ = Ibfxcd ® E^^ ® F**, integration over Ha simply 
causes the substitution hlh^^ ^—> cohw^h^^ , followed by 
restriction of the integration domain from Hp to Hp/HA. 

To do the integral over the latter, we exploit the sim- 
plifications that arise from the limit e ^ 0. What must 
happen is that the integral over the noncompact mani- 
fold Hp/Ha becomes singular as for e ^ 0, thereby 
cancelling the prefactor e^. This singular behavior im- 
plies that the integral is dominated by the contributions 
from the asymptotic region on Hp/HA. We can therefore 
replace the integral by its contraction to the asymptotic 
region. Using the relation h{At — 2TT^)h^^ = At — 27r+ 



we see that 7r_|_ may be replaced by At /2 in the integrand 
with an error that becomes negligible in the limit e — > 0. 
We set Q :— ehAth^^, which satisfies — e^TTf, so 
that = on the contracted (or "light-cone" ) manifold 
that emerges in the limit £ — > 0. The invariant measure 
£ ■ c?/iHp/Ha be shown to contract to the invariant 
measure dQ on the light-cone = 0. Hence, the in- 
tegral giving the single-particle form factor Fp^^ reduces 
to 



/ dQ STrKQ e^'^'^''^'^ 



d 

dt 



dQ e 



-tSTrKQ 



(39) 



To do the integral over Q, we write 



0= E (Q^^) 

i,j=A,B. 



bfxcd 



Eli <E) F* 



and solve the constraint = by putting 



Qaa = ZZ\ Qar = -ZVzTz, 
Qra = ^ZTZ Z\ QflK = -Z^Z, 



with 



Z = E^^ ® 



Zl ^2 
'Z2 Zi 



E 



FB 



cd 



Ci Ci 

C2 C2 



cd 



The invariant measure dQ in these variables is the Eu- 
clidean measure dZdZ^^ = dzi(izi(iz2(iz2i9ci^Ci'^C2^C2 '^'^l" 
tiplied by STr ZZ'^ . Upon making the decomposition 



K 



bfxcd 



Fr 



F' 



i=A,R 



we get STr KQ = STr Z^ {KaZ - ZKr) , and our integral 
becomes 



d 
di 



t=i 



dZdZ^ STr ZZ^ e-tSTr zHKAZ-ZKn)_ 



By the property of perfect grading (i.e. equal number of 
bosonic and fermionic degrees of freedom) , the Euclidean 
measure dZdZ^ is invariant under scale transformations 
Z — > Zj\ft. Therefore, carrying out this transformation 
and then taking the derivative at i = 1, we obtain 



dZdZt STrZZt e-STr^nifA2-ZK,)_ 



This is a Gaussian integral which is easy to calculate. 



The result is the one given in Sect. [V 



We conclude this appendix by sketching the evidence 
in favor of the conjecture (pT|). When proper mean- 
ing is given to the quantities K and Q, the single- 
particle form factor Fp^^ can be written in the form 
(^) for all the cases (A — 1/2,1,2) accessible to our 
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method. The invariant measure dQ takes the general 
form dZ(iZ^(STrZZt)^-i. By doing the same manipu- 
lations as before, we obtain 

= (A - 1) I dZdZt(STrZZt)^-ieST'-^'(^^«-^^^). 

For integral A, the pre-exponential factor {STr ZZ^)^^^ 
can be generated by shifting Kji — > Kn + [k^v /vs) ■ 1 and 
taking A — 1 derivatives w.r.t. w at w = 0. If we switch to 
the notation of Sect. the resulting Gaussian integral 
equals 



/ 



dZdZ'^ g~STr Z^KAZ-ZKn-Zkov/v,) 



in all cases. This is why we believe that (yj) holds true 
in general. 



APPENDIX C: CONSTRUCTION OF 
BOUNDARY TERMS 

We sketch our general procedure for constructing 
the boundary terms that are associated with polar- 
coordinate superintegrals on supersymmetric spaces. 
Consider the simple example of the Euclidean superplane 
E212 with complex supercoordinates (2;, C) and invariant 
Berezin integral 



fdfi={i7Tir' I (dcd-^f{z,Z,U)) 



dzdz. 



The metric of E212, dzdz -I- d^dC, is preserved by trans- 
formations 



if fc e U(l|l). We introduce polar coordinates (fc,r) = 
(j)~'^{z,z,(^X) via the diffcomorphism : U(l|l)/U(l) x 
R+ E212 defined by 



The substitution rule for superintegrals yields 



fd^l 




{f o (l)){k,r)dk r-^dr + nlf] 



U(l|l) 



where we have extended the angular integration to the 
group U(l|l) for convenience. The normalization con- 
stant l/Ani has been absorbed into the Haar-Berezin 
measure dk. 7^[/] is a correction term which is due to 
the radial space R"*" having a boundary. To construct 



it, we exploit the translational invariance of the Berezin 
integral on E212 as follows. We write X := (z,z, C, for 
short and denote the scalar product of two vectors X, X' 
by 

{X,X') = 3? (zz' + CC')- 

If Xq — (zq, zq, Coj Co) is some set of parameters, we define 
a directional derivative -Do/ by 

i?o/ = (^o,grad/). 

From {Dof)dfi — we deduce 



n[Dof] = - f If {Dof 



o 0)(fc, r)dk r ^dr. 



By first decomposing the directional derivative into its 
radial and angular parts, then partially integrating, and 
finally recombining terms, we obtain the relation 

- / iDofo^){k,r)dk = 

(fc-i-Xo,ei>(/o</,)(r, k)dk\ 




U(l|l) 



where ei is the unit vector ei = 



, and k ^ ■ Xq de- 



notes the element Xq rotated by fc ^. We now integrate 
the total radial derivative to get 



n[Daf] = lim - [ (fc-^ • X, 
'■-^o r Ju(i|i) 



o,ei){f o (f)){k,r)dk. 



By Taylor expanding / to linear order in r and doing the 
resulting integral over U(l|l), we obtain the expression 

n[Dof] = (Xo,grad/(0)) = po/)(0), 

from which conclude 

nf] = /(o). 

This elementary result is well-known; for a direct deriva- 
tion see e.g. The above construction has the great 
virtue of using no more than the generic structures at 
hand. It therefore generalizes easily to more complicated 
situations. In particular, it can be used to construct the 
boundary terms associated with the polar-coordinate in- 
tegral The result is as stated in the text. 
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